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The Classical-map Hyper-Netted-Chain (CHNC) technique is a simple method of calculating 
quantum pair-distribution functions, spin-dependent energies, etc., of strongly-interacting uniform 
systems. We present CHNC calculations of charge densities and energies of non-uniform systems, 
viz., quantum dots, and compare with quantum Monte Carlo and density -functional results. Results 
for up to 210 electrons are reported. 
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The Hohenberg-Kohn and Mermin (HKM) theo- 
rems [H of density-functional theory (DFT) assert that 
the one-body density n(f) of an inhomogeneous system 
completely determines its physics. However, DFT uses 
the more laborious Kohn-Sham (K-S) approach Q due 
to the lack of an accurate kinetic-energy functional f^]. 
The Kohn-Sham n{r) of an electron system is: 



The map using the Hyper-Netted-chain (HNC) method 
is named the CHNC. Molecular dynamics (MD) simu- 



n(r) = ^|</..(r)p/.(Wr) 



(1) 



Here are K-S eigenstates with "energies" e^, occu- 
pations factors fi, at the temperature T = 1//3. Here 
V includes spin and other relevant quantum numbers. 
Quantum systems at high temperatures behave classi- 
cally. Then n{r) is given by the Boltzmann form: 



n(f) = nQ{{))exp{-PVKs{r)) 



(2) 



where 7io(0) is a reference density, and Vks is a classical 
Kohn-Sham potential. This suggests that the kinetic- 
energy functional may be side-stepped by (i) the use of 
an equivalent "classical-fluid" at a temperature Tcf for 
the quantum system whose actual physical temperature 
T may even be zero; (ii) use of effective classical pair- 
potentials to mimic quantum effects. In the following we 
discuss such a "classical map" , enabling a great simplifi- 
cation in quantum calculations of interacting electrons. 

All quantum observables are mean values over suitable 
distribution functions, formed by averaging over most 
of the variables in the square of the many-body wave- 
function. The most useful averaged quantities are the 
one-body and two-body distributions. We had already 
demonstrated a simple but accurate classical map for the 
interacting uniform electron fluid (UEF), by presenting 
explicit calculations of spin-polarized pair-distribution 
functions (PDFs) at zero and finite T, for the 3D elec- 
tron liquid 0, H, the 2D-electron fluid [1, S i, i] and 
multi- valley systems [lo| . Fermi-liquid parameters of 
thick-electron layers have also been determined via this 
classical-map technique [ll| . It has been successfully ap- 
plied to hot dense hydrogen and related systems 1^, [l3| . 

The method employs a classical Coulomb fluid whose 
PDFs are determined via classical statistical mechanics. 



lations may also be used where it was called the 



CMMD. The temperature Tq of the "equivalent classical 
Coulomb fluid" is chosen to reproduce the correlation 
energy of the original quantum fluid at T 0. Then 
the PDFs of the classical fluid at T^f = ^J{T'^ + T^) are 
excellent approximations to the PDFs of the quantum 
fluid at T . The so-obtained PDFs are then used in the 
standard way to obtain the energies and other proper- 
ties of the system. Where possible, the accuracy of the 
CHNC results have been demonstrated by comparison 
with quantum Monte Carlo (QMC) or DFT results. 

Here we apply the CHNC to a typical inhomogeneous 
systems, viz., electrons trapped in 2D parabolic poten- 
tials (quantum dots). In standard calculations, if iVe elec- 
trons are in the "external potential" , a suitable basis set 
of Nh functions, with TV^ significantly larger than A^e, 
is selected. A Hartree-Fock (single-determinant) calcula- 
tion is followed by a configuration-interaction (CI) expan- 
sion in Slater determinants. The complexity of the prob- 
lem grows factorially with A^h. It is the electron-electron 
interactions, which make the problem prohibitive. In 
CHNC or CMMD, we treat many-body effects classi- 
cally (i.e., an 0(0) approach), while the non- interacting 
Hamiltonian is treated exactly. Here we summa- 
rize the salient features of the classical-map technique: 
(i) Assignment of a classical-fluid temperature Tc/ to 
the electron system, (ii) Replacement of the Coulomb- 
interaction operator 1/f by a classical "diffraction poten- 
tial" Vcir)= {1 — exp(— fctftr)}/r which accounts for the 
thermal de Broglie length Xjkth of the electron at Tc/. 
(iii) Ensuring that the non-interacting electron PDFs 
with spin polarization C, viz., g^(r, T, ^) are correctly re- 
covered if the Coulomb interaction is switched off. Thus a 
"Pauli exclusion potential" P(r) is used to exactly repro- 
duce the Fermi hole [14] in 5°(r,r, C). We apply CHNC 
to parabolic quantum dots to show that their interact- 
ing inhomogeneous charge densities and energies can be 
readily calculated via CHNC, for arbitrary A"e- 

Quantum dots- Electrons trapped in parabolic poten- 
tials are found in Fermion traps and in the quantum dots 



of nanotechnology[15j. Such quantum dots have been 
studied extensively by several methods 3i 12 |i includ- 
ing QMC [Tl|. Using atomic units {\e\ = h = rrie ^ I), 
the Hamiltonian operator H = Hq + Hint for electrons 
subject to a potential u{r) is given by 



'V2 



Hint — 7: 



^ (3) 



The classical map converts the Hamiltonian operator to 
the classical Hamiltonian Hr. 



1 



2 ^ 



(4) 

Here pi is the momentum of the i-th electron, with spin 
Si. The Coulomb- interaction operator 1/f is replaced by 
the well-known diffraction-corrected classical form Vc{r). 
The operator defining the parabolic confinement in the 
dot is u{r) — ujQr'^/2. It maps to the classical func- 
tion Uc{r), constructed so that the non- interacting den- 
sity n'^{r) is recovered from Uc{r) as a classical distribu- 
tion. The potentials u(f) and Uc{r) differ because the 
quantum system is sensitive to the boundary conditions 
imposed on the eigenstates of Hq 

Thus the essential input to the classical mapping of 
inhomogeneous systems is the non-interacting density 
71" (r). Here we suppress vector notation (unless needed 
for clarity) and consider circular quantum dots. Given 
n'^{r), we seek the classical potential which generated it. 
As this involves the inversion of an inhomogeneous HNC- 
type equation, we proceed indirectly. If the presence of 
each electron did not affect any other electron, the cor- 
responding classical potential Uc{r) is: 



n°(r) = n°(0)exp{-/3ucW} 



(5) 



This equation determines the product (3uc(r), and not 
separately the potential Uc{r), or an inverse tempera- 
ture p. The reference zero of all potentials will be re- 
ferred to the center of the dot. The classical potential 
Uc{x) depends on Ne even though the electrons are non- 
interacting. In effect, they have developed indirect steric 
interactions via the confining potential. 

The confining potential defines a length scale £0 — 
^/ {h / {mujo)) ■ We use the effective mass m* = 0.067 
and the dielectric constant k = 12.4 typical of GaAs. 
These define effective atomic units (a.u.) with a Hartree 
energy of me'^ / {kTiY — 11 meV, and a Bohr radius 
as = K I (me^') ~ 9.79 nm. In Fig. [T]Ja) we show the 
ri°(r) and (3uc{r) for a 20-electron circular quantum dot, 
Ne = 20, with Wo — 3.33 meV [3]. The non- interacting 
density N'^{r) shows shell- filling effects. In Fig. [ija) we 
show a smoothed charge density n° (r) as well. The den- 
sity difference n°(r) — n'^ir) integrates to zero, and may 
be thought of as the result of a steric interaction Ustr {r) 




1 


1 ' 1 ' 1 ' 1 

(b) 




•— • DFT 




CHNC 


- \ 






- n°(r) 




\ a — . n(r)-Pois 




*\ -- CHNC-xc 




\ N =20 




\\ e 






1 


1 1 KenJnJS-L- 



1.6 



1.2 



0.8 - 



0.4 



4 1 



1 1 1 1 1 

N =210 

e 


1 ' 1 
(d) 


N =110 ^\ 




^ N =30 \ 








1 < 1 Vl 


\.^J — 



0.4 



4 6 



10 



FIG. 1; (Color online) Electron densities and potentials for 
a 2D harmonic-oscillator quantum dot with A''e=6,12,20,30, 
and 110 and 210 electrons. The unit of length £0 — 1.89 a.u., 
or 11 nm, with ujo — 0.28 a.u., i.e., 3.3 meV. Panel (a) shows 
the non-interacting n'^{r), and also a smoothed charge density 
ng{r). The corresponding classical potentials Uc{r) and Ucs{r) 
are also shown. The difference potential, Ustr{r), defines a 
"steric -packing" potential which is used to construct a Bridge 
correction B{r) to the HNC equations. Panel (b) shows the 
evolution of n''(r) to the interacting n{r) with the inclusion of 
Poisson, exchange and correlation (xc) terms, to give CHNC- 
xc, which does not contain B{r), and CHNC inclusive of B{r). 
In (c) the final CHNC results (inclusive of bridge terms) for 
A^e=6,12, and 20 are given. In (d) results for A^e=3Q, and 110 
and 220 are shown. The QMC data points were extracted 
from Ref. [13. 



arising from the packing of classical particles into the 
parabolic trap. In the quantum system it arises from 
the boundary conditions on Hq. As many shells are 
filled with increasing N^., and as Coulomb interactions 
come on, this effect becomes unimportant, as seen in the 
interacting-n(r) for the dot with Ne~21Q, in Fig.[T](d). 

The electron-distribution na{r), spin a, in the dot de- 
fines a dot-electron PDF by, e.g., na{r) = rta(0)5da(?')- 
This in turn depends on the electron-electron PDFs gee, 
i.e., explicitely, ga^a' {f,f'). Evaluating this coupled set 
is complex even within classical mechanics. Unlike in 
uniform systems, these PDFs dependent on f,f'. While 
a classical-map-MD approach is the best option, we show 



3 



that simple approximations to the CHNC integral equa- 
tions lead to surprisingly accurate results. 

Our simplified approach is based on constructing 
the density ^^^(r) = na{0)gda{r) while the e-e PDFs 
ga,a'{r,r') are evaluated from an average-density ap- 
proximation. This saves us from solving a set of coupled 
HNC equations. The approximations proposed are: (i) 
replacing the e-e PDFs ga,a' (r, r") by the PDFs of a uni- 
form slab (USB) of average density, n, (ii) determining 
the USB density, viz., n from < n{r)n(r) >/< n{r) >, as 
in Refs. (iii) using n and the equations of Ref. 0| 

to determine Tq and other UEF parameters needed for 
the CHNC, (iv) calculating the charge density ^^^(r) in 
the quantum dot via a simplified HNC-like classical in- 
tegral equation inclusive of a bridge term B^^ir). For 
brevity of presentation, we use a spin-unpolarized sys- 
tem [C, — 0), and write gde{r) for the PDF defining the 
charge density in the dot. The basic CHNC equation for 
the electron density in the inhomogeneous system is: 



n(r) = n(0) exp [-f3uc{r) + V„i-b] 



(6) 



This is a Boltzmann distribution for the pair-potential 
Uc{r) inclusive of its many-body correction Vm-t- This 
includes mean- field and correlation effects. The classi- 
cal mean-field term is just a Poisson potential. We also 
need the correlation potentials beyond mean-field from 
classical mechanics. Hence we rewrite Vm-b in terms of 
the Nodal function N{r) and the bridge function B{r) of 
HNC theory, remembering that "exchange" already ap- 
pears in the classical map as an effective pair-potential 
("Pauli potential") between like-spins. The "bridge" dia- 
grams brings in irreducible (three-body and higher) pack- 
ing effects beyond simple HNC. 



N{r)+Bde{r) 
P[Vp{r) + V^{r) + V,{r)] 

n{r')dr' /\r — r'\ 



(7) 

(8) 

(9) 
(10) 

14(r) = I n{r')d? \og{gee{r.r')} - K,{r,?)][ll) 



Vm-b 

N{r) 
Vp(r) 

V;(r) = / n{r')df'P{f,r') 



The nodal term N{r) has been decomposed into a Pois- 
son potential Vp{r), an "exchange potential" Vx{r) aris- 
ing from the Pauli exclusion potential P{f,r'), and a 
classical correlation potential Vc{r). This is a standard 
analysis based on the Ornstein-Zernike equation, and was 
already discussed in Ref. 2]J. The expression for Vc{r) 
is Eq. (3.4) given there, and hee{r,r') is the total corre- 
lation function. Vp has been written with the Coulomb 
potential 1 /r rather than with the diffraction correction 
(DC), since we found that the numerical effects of DC 
are negligible for the gde{r) of the quantum dots studied 
here (n.b. the DC is needed in the CHNC evaluation of 
the geei^) at the uniform-slab density n). 



We approximate gee{r,r') by gee{\r — r'\) of the uni- 
form slab at density n. The n in our quantum dots 
are found to have Wigner-Seitz radii ~ 1 and hence 
the Coulomb correlation potential Vc(r) is small, while 
the Pauli exclusion effect is large. Here P(r, r ') is re- 
placed by P{\r — r'\) at the uniform density n. This 
is a universal function oi x = r/vs- It can be fitted 
to {— a21og(2:) -|- 6}/(l + ex), or to the simpler form 
ai/(l + a2x), with ai=10.1186, 03 = 3.69352, for 2D 
systems at T = 0. Replacing geeir^f') by its uniform- 
density value leads to the question of the appropriate 
form for n(r ') which now has to play the role of a n{r, r'). 
The usual simple choices (e.g., n(r ') {n{r) + n{r' j)/2, 
or even n{r) give similar results, but the replacement 
n(f') — > n, i.e., completely by a uniform slab, is too 
drastic. 

The results obtained from the self-consistent solution 
of Eq. [6] are shown in Fig. [T] Panel (a) shows the in- 
put potential (3uc{r) based on the non-interacting den- 
sity nP{r). This n'^(r) gets modified by the interactions. 
The average density n, the corresponding inverse temper- 
ature (3, potentials Vp, Vxc etc., were calculated from n(r) 
in each iteration, with the total number of electrons fixed 
to Nq. In panel (b) the self-consistent n(r) obtained as we 
successively add Poisson (curve with boxes), exchange- 
correlation (dashed red curve), and Bridge corrections 
(solid blue curve) is found to converge to the benchmark 
results from QMC and DFT (solid green circles). Panels 
(c), (d) shows the full CHNC results and QMC and/or 
DFT results for 7Ve=6,12,20, 30. These have been done 
using the quantum temperature Tq defined by Eq. 5 of 
Ref. ^ , which assignes a ^ to a given . For high electron 
densities when Ts < 1, the Buluty-Tanatar(BT) mapQ 
seems to be more accurate. Calculations for 7Ve=110, 210 
etc, using a BT-type map are given in panel (d), but no 
published QMC or DFT results are available.. 

An oscillatory structure in geeij") occurs even in uni- 
form fluids, when interactions are important, and is 
well understood. The geeif) of the uniform fluid at 
strong coupling could be accurately recovered on includ- 
ing bridge contributions Bee{r) to the HNC, as shown 
in Ref. [a, [3] ■ There the particle-packing theory of the 
hard-sphere fluid could be used, since the PDF is not 
too sensitive to the details of the bridge interaction. In 
the classical map of the quantum dot, packing effects are 
dominated by the steric crowding effect of the confining 
potential. This steric potential Ustr is already available 
to us in the charge distribution ri°(r) of Ho. We assume 
that the difference between the smoothed distribution 
n^{r), and n'~'{r), Fig[Tl^a), corresponds to the effect of 
Ustr{f)- Then, converting a charge distribution into a 
potential via the classical map, we write 

pUstrir) = -log{n°(r)/n°(0)}-logK(r)/nO(0)} 



(12) 



The equation for Bde{r) reflects the rescaling of the 
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TABLE I: The Exchange- Correlation and kinetic energies 
evaluated from the densities (Fig. 1), i.e., CHNC n(r), and 
from the DPT n(r) of Ref. Tt^. A^e is the number of electrons. 
The energy unit is luq. 



ing to a classical molecular-dynamics implementation of 
the method [25!. The author thanks Mario Gattobigio for 
providing the DFT n(r) and energies of Ref [l3|- 
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uniform-slab density "nP to n due to interactions, chang- 
ing the scales of the parameters /S", r°, etc., to (3, 
etc., of the final self-consistent density, thus rescaling 
the steric potential of the non-interacting system. The 
numerical factor 7 is set to 1.5. This simple model 
avoids the complex microscopic calculation of a bridge 
correction, and is seen to be justified a postiori. It does 
not appeal to any parameterizations outside the prob- 
lem. In fig. (He) we show comparisons of the CHNC n(r) 
for Ne =6,12,20, with QMC results. DFT results for 
Ne — 20 are shown in fig. [Hb). Panel (d) shows the 
CHNC resuhs for iVe=30, 110 and 220 electrons. In the 
last two cases we do not have microscopic calculations for 
comparison with CHNC. The calculation of the interact- 
ing density n(r) for arbitrary N^, at finite temperatures, 
finite values of ( or finite magnetic fields pose no addi- 
tional difficulty in CHNC. 

The charge distributions of CHNC can now be used for 
the total energy E, which involves the confinement en- 
ergy Ec, the Possion energy Epo^, E^c, and i?fej„. The 
exchange-correlation and kinetic contributions are the 
quantum mechanically sensitive, "difficult" terms. The 
simplest approach is found to be adequate for 2D quan- 
tum dots. That is, we use the LDA (local-density ap- 
proximation), with the known 2D exchange-correlation 
energy functionals [22] . The success of the LDA for the 
2D kinetic energy has been noted by van Zyl et al. 2J], 
and also by Koivisito et al 2]. A comparison of our xc- 
energies with those from DFT are given in Table I. 

In conclusion, we have presented classical-map calcu- 
lations for a 2D inhomogeneous system of interacting 
electrons, viz., quantum dots, which are in good agree- 
ment with microscopic calculations where available. This 
method requires no basis sets, no evaluation of matrix 
elements etc. It is an order-zero, viz., 0(0) approach 
independent of the number of electrons. Similar applica- 
tions to atomic systems, ("naturally occurring quantum 
dots"), are clearly feasible. The uniform-slab approxima- 
tion, and the need to model the bridge-diagram correc- 
tions to the CHNC equations may be avoided by resort- 
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